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Abstract 

This paper investigates the problem of Gaussian approximation for the wireless multi-access 
interference distribution in large spatial wireless networks. First, a principled methodology is 
presented to establish rates of convergence of the multi-access interference distribution to a 
Gaussian distribution for general bounded and power-law decaying path-loss functions. The model 
is general enough to also include various random wireless channel dynamics such as fading 
and shadowing arising from multipath propagation and obstacles existing in the communication 
environment. It is shown that the wireless multi-access interference distribution converges to the 
Gaussian distribution with the same mean and variance at a rate where A > is a parameter 
controlling the intensity of the planar (possibly non-stationary) Poisson point process generating 
node locations. An explicit expression for the scaling coefficient is obtained as a function of 
fading statistics and the path-loss function. Second, an extensive numerical and simulation study 
is performed to illustrate the accuracy of the derived Gaussian approximation bounds. A good 
statistical fit between the interference distribution and its Gaussian approximation is observed 
for moderate to high values of A. Finally, applications of these approximation results to upper 
and lower bound the outage capacity and ergodic sum capacity for spatial wireless networks 
are illustrated. The derived performance bounds on these capacity metrics track the network 
performance within one nats per second per hertz. 
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I. Introduction 

A. Background and Motivation 

Wireless communication technologies have been evolved at an unprecedented pace over the 
last decade. This has, in turn, given rise to the birth of many next generation wireless systems 
including 4G networks, femtocells and cognitive radio networks. The main design philosophy 
underpinning most of such emerging classes of wireless systems, if not all, is the more proficient 
utilization of wireless spectrum than ever before to accommodate high volume of data traffic from 
increasingly overcrowded user populations within a given frequency band as well as to support 
increasingly more data-rate-intense multimedia applications over the wireless links. Interference 
mitigation and characterization are the primary design challenges to overcome for achieving this 
goal and for meeting the high target data rates (e.g., 1 Gbit/s for low mobility and 100 Mbit/s 
for high mobility environments for 4G networks [JT]) set for future wireless networks consisting 
of many interfering links. 

However, even as a first step to characterize network performance and relevant performance 
metrics, computation of the exact wireless multi-access interference (WMAI) distributions turns 
out to be a mathematically intractable problem in most practical scenarios. This is mainly because 
the level of WMAI at a receiver node in a wireless network is a random function of the overall 
network geometry as well as numerous other wireless channel dynamics such as path-loss, fading, 
shadowing and so on, e.g., see ||2]-||4). This motivates the search for tight bounds on the WMAI 
distributions that can accurately track the interference behavior [0-|El. 

This paper focuses on the statistical characterization of WMAI for spatial wireless networks, 
and establishes tight Gaussian approximation bounds for the WMAI distributions. To this end, the 
underlying spatial stochastic process determining transmitter locations is assumed to be Poisson 
but not necessarily stationary. The signal power attenuation due to path-loss is modeled by 
means of a general bounded and power-law decaying path-loss function. Other random wireless 
channel dynamics such as fading affecting the received signal power are also accounted for in 
the employed signal propagation model. 

Our main contribution is the derivation of the rate of convergence of the WMAI distributions to 
the Gaussian distribution with the same mean and variance, which is formally stated in Theorem 
[TJ Briefly, this rate is equal to where A is a modeling parameter enabling us to control 
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the "intensity" of the planar Poisson point process (PPP) generating transmitter locations (i.e., 
see Section [TT] for details), and c:l4 R + is a positive function which depends on the shape 
of the path-loss function and the point x G R at which we want to estimate the value of the 
(centered and normalized) WMAI distribution function. c(x) approaches zero for large absolute 
values of x at a rate \x\~ z . This behavior makes the derived bounds on the tails of the WMAI 
distributions tight for any given value of A. Moreover, the supremum of c(x) over x is a small 
constant, which allows us to obtain uniform rates of convergence as a function of A. 

An extensive simulation and numerical study is performed to illustrate the accuracy of derived 
Gaussian approximation results estimating WMAI distributions for specific instances of both 
stationary (i.e., see Section HVT) and non- stationary (i.e., see Appendix iDl) PPPs. It is observed 
that the Gaussian distribution with the same mean and variance can accurately track the be- 
havior of the WMAI distribution even for small values of A. As predicted by our bounds, the 
approximation accuracy increases further as A increases, and the WMAI distribution becomes 
almost undistinguishable from its Gaussian approximation when A is around 10. On the other 
hand, random fading effects have an adverse impact on this Gaussian behavior, and the bounds 
given in Theorem Q] become looser when the fading distributions have larger dynamic ranges. 

The utility of the derived Gaussian approximation bounds on the WMAI distributions is also 
illustrated by characterizing wireless network performance in tangible communication scenarios. 
In particular, tight upper and lower bounds on the outage capacity for a given victim link in 
a spatial wireless network and those on the sum capacity for spatial wireless multiple-access 
networks are obtained. The approximation accuracy of the derived performance bounds on these 
network capacity measures lies within one nats per second per hertz for moderate to high values 
of A. For high values of A, our bounds become very tight, and they almost coincide with the 
outage capacity (i.e., see Subsection IV- Al) and the sum capacity (i.e., see Subsection IV-BI) . 

B. Related Work and Paper Organization 

Wireless networks are often interference limited due to the broadcast nature of the wireless 
medium. Hence, the statistical characterization of interference for large spatial wireless networks 
has been a key area of research for more than a decade. Many of the latest developments in the 
field are summarized in the recent review articles [0-0]. 

Historically, the early efforts for characterizing the structure of WMAI in wireless networks by 
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using stochastic geometry can be traced back to as early as 1978 [8|. Sousa et al. applied similar 
techniques, e.g., [9] and [10], in the 1990s to assess the performance of spread spectrum wireless 
networks as well as to find optimum transmission ranges in these networks. Subsequently, several 
approximation techniques appeared in the field to approximate the WMAI distributions, and 
then to use these approximation results for the network performance analysis, e.g., Q, [fTTI and 
|fT2l. However, except for one special case Q, i.e., the case in which the signal power decays 
according to the unbounded power-law decaying path-loss function t~ 4 , there are still no closed 
form expressions available for the WMAI distributions. 

Therefore, it becomes necessary to resort to numerical methods to calculate the WMAI 
distributions by modeling WMAI as a power-law shot-noise process and then inverting its 
computed characteristic function. Even though algorithmic perspectives based on fast Fourier 
transforms lfl3l to numerically compute power-law shot-noise distributions and densities are 
promising, they are of limited interest and importance in the context of wireless networking. The 
numerical computation cannot provide closed form expressions revealing structural dependencies 
between the WMAI distributions and network design parameters to assess the wireless network 
performance under candidate/existing wireless communication technologies. 

Alternative approaches in the field include various approximation methods based on geomet- 
rical considerations and distribution approximation techniques ED-ATI, LePage series [fTTI or 
Edgeworth expansion [fT2j to estimate the WMAI distributions. The main objective is to obtain 
simple but insightful upper and lower bounds on the related performance quantities of interest 
in spatial wireless networks by utilizing such approximations. In particular, our results in this 
paper are close to those in [|2]-[]7j in that the authors of these works have also obtained upper 
and lower bounds on the normalized WMAI distributions with provably small gaps between the 
computed bounds and the actual WMAI distribution!] 

In contrast to [5] as well as most other earlier works in the field such as [9] and [10], one 
distinctive aspect of our analysis in this paper is that we work with general bounded and power- 
law decaying path-loss functions to calculate the WMAI distributions. The underlying motivation 
for working with bounded path-loss models is recent findings pointing out that the unrealistic 

'lb be more precise, the authors in J5] focused on the distribution of the normalized inverted signal-to-interference-ratio, 
which is, in essence, the same as computing the WMAI distribution. 
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singularity of the unbounded path-loss model at leads to unexpected deviations on the final 
computed WMAI distributions, e.g., see lfl4l and [fT5l . When compared with the results reported 
in (6]| and Q, we derive tighter Gaussian approximation bounds for more general spatial node 
distributions including non-stationary PPPs and for more general signal propagation models 
including fading and shadowing effects. 

Since interference in a wireless network is a specific instance of a shot-noise process, the 
results of this paper are also related to a more general body of work on shot-noise processes 
such as lfT6ll and ifTTll . The paper lfT6l established many properties for power-law shot-noise 
processes on the line such as its moment generating functions, moments and cumulants. For a 
very specific bounded and power-law decaying impulse response function driving the power-law 
shot-noise process, the convergence of the amplitude distribution of the power-law shot-noise 
process to the Gaussian form was briefly mentioned in [fl6l . but without any formal proof for 
this convergence and without establishing rates of convergence. 

Our results are closely related to those reported in [17] since the authors in [fTTl also established 
the uniform rates of convergence for the amplitude distribution of the shot-noise process. When 
compared with [fTTl . our results are different than those of [fT71 in three aspects. Firstly, we 
give both uniform and non-uniform Berry-Esseen types of bounds on the WMAI distributions. 
The non-uniform bounds allow us to tightly approximate the WMAI distributions for extreme 
interference values even for small finite values of A. Secondly, error terms appearing in our 
uniform Gaussian approximation formulas are eleven times smaller than those in ifTTI . This 
implies much tighter performance bounds for characterizing various performance measures for 
spatial wireless networks. Finally, this paper introduces a principled and simpler methodology 
to establish Gaussian approximation results. Hence, the distribution approximation techniques 
presented here are expected to find other potential applications in the analysis of emerging 
wireless network architectures. 

The remainder of the paper is organized as follows. Section U describes the network configu- 
ration along with our modeling assumptions. In Section [Till we establish the main Gaussian 
approximation result for the WMAI distributions. Section |IV] illustrates the applications of 
the derived Gaussian approximation bounds for stationary PPPs, and presents an extensive 
numerical and simulation study to verify these bounds. Similar analysis is also performed for 
a non-stationary PPP in Appendix |Dj Section |V] presents further applications of our Gaussian 
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approximation bounds to characterize outage capacity and sum capacity for spatial wireless 
networks. Finally, Section |VT] concludes the paper. Most of our proofs are relegated to appendices 
for the sake of fluency of the paper. 

II. Network Model 

We consider a planar network in which transmitters are distributed according to a planar PPP 
with mean measure A (alternatively called: intensity measure), denoted by $a, over M 2 . Here, 
for any (Borel) subset A of M 2 , A (^4) gives us the average number of transmitters lying in 
A. We will assume that A (A) is locally finite, i.e., A (A) < oo for bounded subsets A of IR 2 , 
and A (M 2 ) = oo, i.e., there is an infinite population of transmitters scattered all around in M 2 . 
The location of the kth transmitter is represented by X k . We will often represent the transmitter 
location process $ A as a discrete sum of Dirac measures as $ A = J2k>i &x k -, where 8 Xk (A) = 1 
if Xk G A C M 2 , and zero otherwise. -XV s can be interpreted as points (or, atoms) of $a, and 
therefore, we also use the notation Xk G $a to symbolize this interpretation. For a point x G M 2 , 
i = 1,2, represents the zth component of x. B(x,r) represents the (planar) ball centered 
at x with radius r. 

We consider the case in which all transmitters transmit with the same power P, i.e., the case 
of a non-power controlled wireless network. For the signal power attenuation in the wireless 
medium, we consider a bounded monotone non-increasing path-loss function G : [0, oo) 
[0, oo), which asymptotically decays to zero at least as fast as t~ a for some path-loss exponent 
a > 2. In addition to the signal power attenuation due to path-loss, it is assumed that transmitted 
signals are also corrupted by fading. The random (power) fading coefficient at transmitter k is 
given by H k , and is assumed to be independent of ^aJ 2 ! The received signal power at a distance 
is given by PH k G (tk) for transmitter k. H k 's are independent and identically distributed (i.i.d.) 
with a common density q(h), h > 0, and have finite first, second and third order moments, 
i.e., J °° h k q(h)dh < oo, k = 1,2 and 3. The first, second and third order moments of the 
fading coefficients are denoted by m H , m H 2 and m H i, respectively. We note that the employed 

2 For simplicity, we only assign a single fading coefficient to each transmitter. In reality, it is expected that the channels 
between a transmitter and all potential receivers (intended or unintended) experience different and possibly independent fading 
processes. Our simplified notation does not cause any ambiguity here since we focus on the total interference power at a given 
arbitrary position in E 2 in the remainder of the paper. 
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signal propagation model is general enough that H k 's could also be thought to incorporate 
shadow fading effects due to blocking of signals by large obstacles existing in the communication 
environment although we do not model such random factors explicitly and separately in this 
paper. We also note that our model is general enough to include the case where transmitted 
signal powers are also random due to power control, i.e., we just need to scale the fading 
process with transmitted signal powers in this case. As it is common in most of the earlier 
works, we will assume a contention based medium access control (MAC) layer such as ALOHA 
mediating node transmissions, and giving rise to the observed distribution of active transmitters 
over R 2 . 



We place a test receiver node at an arbitrary point X Q = \^Xo ,Xo j G M 2 , and consider 
signals coming from all other transmitters, whose locations are given by $a, as interference to 
this test receiver node. Without loss of generality, we focus on the distribution of WMAI seen 
by such a test receiver node placed at arbitrary X D for the rest of the paper. If $a is stationary, 
the interference statistics seen from any other point in M 2 will be the same. 

The level of WMAI at X Q depends on the distances between the points of <3>a and X Q . Hence, 
the transformed process J2k>i $T(x k )> whose points lie on the positive real line, is of particular 
importance to derive interference statistics at X Q , where T : M 2 i— > R is given by T (x) = 
\\x — X \\ 2 = J — Xo 1 ^ + — Xo 2 ^ . The mean measure of the transformed process 
is equal to AoT -1 , where T _1 (A) = {x G K 2 : T (x) G A} for all iCE, and AoT -1 (A) has 
the same interpretation above, i.e., it gives us the average number of points of the transformed 
process in A (see the Mapping Theorem on page 18 in ED). We assume that A o T _1 has a 
density in the form A o T _1 (A) = A f A p(t)dt such that p(t) = O (t"" 1 ^) as t -)> oo for some 
e > 0. This assumption on the rate of growth of p(t) is necessary to ensure the finiteness of 
WMAI at X Q . Here, A is a modeling parameter, which can be interpreted as the transmitter 
intensity parameter, that will enable us to control the average number of transmitters lying in 
A and interfering with the signal reception at the test receiver at X a . 

The level of WMAI at X Q is equal to I x = J2k>i PH kG (T (X k )) . I x is a random variable 
since transmitter locations and associated fading coefficients, i.e., {{Xk, Hk) , k > 1}, are random 
variables. Therefore, different node configurations and fading states result in different levels of 
interference at the test receiver. In the next section, we will show that the distribution of I\ can 
be approximated by a Gaussian distribution. 
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A Note about Notation: We use boldface and calligraphic letters to denote vector quantities and 
sets, respectively. | ■ | notation is used to measure the magnitudes of scalar quantities, whereas 
|| • || 2 notation is used to measure the Euclidean norms of vector quantities. As is standard, 
when we write f(t) = O (<?(£)), f(t) = (g(t)) and f(t) = o (g(t)) as t — > t for two positive 
functions f(t) and g(t), we mean limsup^ ^ < oo, liminf^ g| > and lim t _^ ^| = 0, 
respectively. f(t) is said to be (g(t)) as t — > t if fit) = O (g(t)) and fit) = Q (g(t)) as 
t — > t . With a slight abuse of notation, we sometimes use I\(P) notation to represent the 
level of WMAI at X Q when we need to put emphasis on transmission powers to explain some 
performance results in Section [V] 

III. WMAI Distribution and Rates of Convergence to the Gaussian 

Distribution 

This section presents calculations for approximating the WMAI distributions as a Gaussian 
distribution, and establishes the rates of convergence for this approximation as A grows large. 
By using Laplace functionals of Poisson processes (see [fT8l for details), we have the following 
Laplace transform for I x : 

(POO POO \ 
-Xj J (1 - e - sPhG{t) ) q(h)p(t)dtdhj , 

where s > 0. The next lemma shows that I\ has a non-degenerate distribution. 

Lemma 1: For all s > 0, J °° J °° (l - e - sPhG ^) q(h)p(t)dhdt < oo. 

Proof: Please see Appendix [A] ■ 

We will need some auxiliary results to prove the main approximation result of the paper. The 
proofs of these auxiliary results are relegated to the appendices at the end for the sake of fluency 
of the paper. The next lemma shows that the distribution of I\ can be approximated as a limit 
distribution of a sequence of random variables /„, i.e., I„ A I\ as n — > oo. 

Lemma 2: For each n, let Ui n , . . . , Ur& n i n be a sequence of i.i.d. random variables with a 
common density f(t) = ^x^-l{o<t<n}, where A n = A J™ p(t)dt and [•] is the smallest integer 
greater than or equal to its argument. Let 

rA»i 

I n = Y,PH k G{U k)n ). (1) 

k=l 

Then, I n converges in distribution to I\, which is shown as I n — > I\, as n — Y oo. 



s 



Proof: Please see Appendix [Bj ■ 
The next lemma shows that the mean value and variance of I\ can also be approximated by 
the mean value and variance of /„. 

Lemma 3: Let I n be defined as in (OQ). Then, 



and 



lim E [I n ] = E [I x 



lim Var [I n ] = Var [I x 



Proof: Please see Appendix |Cj ■ 
Lemma 4: Let £1, . . . , £j be a sequence of independent and real- valued random variables such 
that each of which has zero mean and Y%=i E l^f] = 1- Let % = Yjl=i E [l&f] ■ Then, 



Pr 



8=1 



< x mm 0.4785 



31.935 
1 + Ixl 3 



for all x e M. 

Proof: The main theorem stated in page 1 in [19] implies that 
X0.4785. Theorem 1 in [|20l after some simplifications implies that 



Pr< 




1 £i ^ X 


> - 


Pr< 




=1 & < Z 





< 
< 



X 



31.935 

T+RF 



. Combining these two bounds completes the proof. 



By using the above assisting lemmas, the main approximation result of the paper is established 
in the next theorem. 



Theorem 1: For all x 6 



Pr 



Eft 



'Varft 



< x 



< 



c(x) 

71 



(2) 



t 

where ^f(x) = -k= f x e~^~dt, which is the standard normal cumulative distribution function 



(CDF), and c(x) 

[m H 2) 

Proof: We let £ fc>n = 
are as defined above, 



f™G 3 (t)p(t)dt 



(/0 00 &(t)p(t)dty- 



min [ 0.4785, ™ 



PH k G{U k , 



for n > 1 and 1 < A; < |"A n ~|, where Uk, n and A n 
E [PH k G (U k>n )} and <7 n = y/\lar\I n }. We note that E [f fcjf J = 
and Efc=i' E = 1- Hence, the collection of random variables {C&,n}[=^ is m the correct 

form to apply Lemma HI We need to calculate Xn = E[=™^ E [|£fc, n | 3 ] to finish the proof. We 
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can upper bound this summation as follows. 



Xr, 



i 1 ■-< 

° n k=i 

® n 



E [P 3 H 3 G 3 (U 1>n ) + 3P 2 H 2 G 2 (U hn ) mi, n + 3PH.G (U hn ) m\ n + m 3 J 



3 [A w ] Pm H m 2 n X f n 

A n<T 3 n 



G 3 (t)p(t)dt + 



3 [A n ] P 2 m H 2m hn X f n 



G(t)p(t)dt + 



Ko 3 n 



G 2 (t)p(t)dt 



Observing m lri = o(l) and m\ n |~A n ] = o(l) as n — y oo, i.e., see Appendix and the 
convergence of a 2 to Var as n — y oo, z'.e., see Lemma [31 we have 

P 3 ra#3A 



lim sup Xn < 



G 3 (t)p(t)dt. 



(Var[/ A ])Wo 

Substituting the expression given for Var [J A ] in (fT9b in Appendix O we can rewrite © as 

mm G 3 {t)p{t)dt 



lim sup Xn < 



V^m^)* (J °°G 2 (t)p(i)dt) 



(3) 



(4) 



By using Theorem |4l we have 



k=l 



< Xnmin 0.4785, 



31.935 
1 + Id 3 



(5) 



for all n > 1 and x G E. Lemmas |2] and [3] imply that 



E f d -<A 

?fc,n — > ~ 



E[/a 



as n — y oo. 



v'VaTTT 

Hence, using © and taking the lim sup of both sides in ©, we have 



lim sup 



Pr<J ^Zk, n <xj-*(x) 
h - E [/a] 



fe=i 



Pr 



< x 



< 



1 mm J °° G 3 (t)p(t)dt 



m H .y (j™ G\t)p{t)dt) 



min 0.4785, 



31.935 
1 + Ixl 3 
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which completes the proof. ■ 
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Firstly, we note that m H 3 > (m H 2) 2 by Jensen's inequality, with equality for deterministic 
fading coefficients. Therefore, the bounds given in Theorem \T\ become tighter if the fading 
distribution is more concentrated around a point, i.e., fading distributions with more restricted 
dynamic ranges. Secondly, the Gaussian approximation bound derived in Theorem \T\ is a com- 
bination of two different types of Berry-Esseen bounds, one of which is a uniform bound and 
the other one is a non-uniform bound. The non-uniform bound is designed to be tight for large 
values of \x\. On the other hand, the uniform bound is tighter for moderate values of \x\. This 
point will be further illustrated in detail in the next section by using tangible examples for the 
transmitter location process $a- It should also be noted that these convergence rates depend on 
the reference point X at which we measure the interference power level. One easy corollary 
of Theorem \T\ is the following uniform approximation bound. 

Corollary 1: 



sup 



pr , h E[J A ] < , _ 



* 7T (6) 



where c = 0.4785 m " 3 3 fJ^M)^ 



m H^y do 00 &{t) P {t)Aty 



IV. Gaussian Approximation Bounds for Stationary PPPs 

In this part, we take $a to be a stationary PPP on M 2 with intensity A [nodes per unit area], 
i.e., A (A) = \v (A) and v (A) is the area of the set A, to numerically illustrate the validity of 
our bounds derived in Section [Till Further applications to approximate the WMAI distribution 
for a non-stationary PPP are demonstrated in Appendix |D] 



A. Theoretical Results 

There are several equivalent ways to represent $a in this case. The most convenient repre- 
sentation for our purposes in this paper is the one obtained by transforming and marking (see 
[fT8l for the details of marking and transforming of PPPs) a stationary PPP with intensity 1 on 
[0, oo), which is given as 

$ A A X^5, rr- frT \, (7) 
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where X = Y means two random variables X and Y are equal in distribution, U k s are i.i.d. 
random variables uniformly distributed over [0,27r], and T k = Yli=i^ where Ei's are i.i.d. 
random variables with unit exponential distribution. This representation allows us to take X k = 

(\f^ cos ' sm ' an( * 11^- ^11 2 = \/&- ^ e assume mat me test receiver is placed 
at the origin, i.e., X a = 0. By using Poisson process transformations one more time, one can 
further show that the distances between the origin and the points of <£>a form a PPP on [0, oo) 
with mean measure A o T _1 ([0, t]) = \%t 2 , and the density p\(t) = 2\irtl^ t > y Hence, the 
WMAI distribution for stationary PPPs can be approximated as in the following theorem. 

Theorem 2: Assume $a is a stationary PPP with intensity A transmitters per unit area. Then, 
for all x G R, 



p r <! h , e[/a] <x \ - *(./■; 
Var[/ A ] " 



c(x) 

where = -4- P e"#dt and c(x) = -4- ^ JTgfffl^ min f _ 47g5 3035 
Proof: Directly follows from Theorem \T\ after substituting 27r£l{ t >o} forp(t). 



5. Simulation and Numerical Study 

Now, we present our numerically computed Gaussian approximation bounds and simulation 
results confirming the theoretical predictions above. We will use two different path-loss models 
Gi(t) = nj^w and G 2 (t) = with various values of a. Similar conclusions continue to hold 
for other path-loss models. 

In Fig. [Q we plot our numerically computed bounds for the WMAI distributions both with 
(bottom figures) and without (top figures) fading. Nakagami-m fading model [|2T| with unit mean 
power gain and m parameter set to 5 is used to model the likely fading effects existing in the 
communication environment. We set the path-loss exponent a to 4 in Fig. [Q 

We observe two different regimes in our computed bounds for the WMAI distributions in 
Fig. ED For the moderate values of the centered and normalized WMAI, i.e., V e[Ja1 , our 
uniform Berry-Esseen bound °- 4785 ^° ° ^ tdt s gives better upper and lower bounds around 
the normal CDF for the interference distribution. On the other hand, for large (i.e., greater 



than 4) absolute values of the centered and normalized WMAI, our non-uniform Berry-Esseen 

31.935 f™G3(t)tdt 

^(i+M 3 ) G 2 (t)tdt y 



bound /^;? 3 ^ — — ^ — — 3- becomes a better estimator for the interference distribution. Our 
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Fig. 1. Upper and lower bounds on the centered and normalized WMAI CDFs for the path-loss functions Gi(t) = (i+t)° 
(lefthand side figures) and G2{t) = (righthand side figures). The effect of fading is also illustrated in the bottom figures 

by assuming Nakagami-ra fading with m parameter set to 5. (a — 4) 



bounds can be used to bound the probability of outage, outage capacity and etc. in a wireless 
communications setting as illustrated in Section |V] We remark that ([8]) is given with an almost 
eleven times larger constant (i.e., 2.21 rather than °i^ 5 ) in fTTvTl (i.e., see Equation 7.1 in ifTTl '). 
Hence, such bounds will be at least eleven times sharper than those based on the convergence 
results in [fTTIl . 

For any fixed value of A, the gap between the upper and lower bounds vanish at a rate 
O (\x\~ 3 ) as x — > oo, i.e., as the interference power increases. When A increases, the upper and 
lower bounds approach the normal CDF at a rate 4=, and we start to approximate the WMAI 
distribution as a Gaussian distribution increasingly more accurately. When the upper and lower 
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bounds on the WMAI distribution are compared for different path-loss models, we see that they 
become tighter for G^Ct). This is because the path-loss dependent constant, i.e., /o G ^ tdt g ; 
appearing in Theorem |2]is smaller for than that for G\(t) {i.e., see Table U below). 

For the simulation study, we built a C-Simulator in order to perform Monte-Carlo simulations 
deriving WMAI CDFs numerically. In Fig. [2l the simulated WMAI distributions are plotted for 
Gi(t) and G%{t) for two different choices of a, i.e., 3 and 5, under various values of A. Fading 
effects are also illustrated in the bottom figures by using the Nakagami-m fading model with 
unit mean power gain, m parameter is set to 1, which is the Rayleigh fading case. Similar results 
continue to hold for other bounded path-loss models and different values of a greater than 2. For 
each different choice of the path-loss model, path-loss exponent and node intensity, we compute 
the interference power level at the origin for 10 4 random node configurations in order to estimate 
the WMAI distributions. 

We focus on small to moderate values of A to illustrate the Gaussian convergence result 
predicted by Theorem El and to understand the effect of small values of A on the WMAI 
distributions. As observed in Fig. |2] the deviations between the normal distribution and the 
simulated WMAI distributions are prominent for sparse networks, i.e., A = 0.1. On the other 
hand, the match between the normal CDF and the simulated WMAI distributions is promising 
for small to moderate values of A, i.e., A = 1 and A = 10. In the non-fading case, for example, 
the match between the simulated distributions and the normal distribution is almost perfect for 
both path-loss models and path-loss exponents when A is around 10 nodes per unit area. Even 
when A is around 1, it is still very good. These observations in conjunction with Theorem |2] 
illustrate the utility of the Gaussian approximation of the WMAI distributions for small and 
large values of A. They also indicate the potential to further tighten the upper and lower bounds 
that we derive in this paper. 

For communication environments with fading, the deviations between the simulated distribu- 
tions and the normal distribution tend to be larger when compared to those for communication 
environments without fading. This is also in accordance with our bounds in Theorem |2] That 

3 

is, the factor m H z/ (m H 2)^ appearing in Theorem |2] is always greater than 1 due to Jensen's 
inequality, as mentioned above, which implies looser approximation bounds in the presence 
of random fading effects. On the other hand, as m increases, the fading process corrupting 
transmitted signals becomes more deterministic, and therefore our bounds become tighter, and 
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2. Comparison of the simulated centered and normalized WMAI CDFs with the normal CDF for the path-loss functions 

(t) ~ (i+t) a ( ie fth an d side figures) and G2{t) = j^jw (righthand side figures). The effect of fading is also illustrated in 

bottom figures by assuming Nakagami-m fading with m parameter set to 1. 



we start to observe better matches between the simulated WMAI distributions and the normal 
CDF. 

When the effect of small and large values of a on the WMAI distributions is analyzed, it is 
seen that the match between the simulated WMAI distributions and the normal distribution is 
slightly better for small values of a (e.g., for A = 0.1 and G2(t) without fading in Fig. |2l the 
maximum deviation between the simulated WMAI distribution and the normal distribution is 0.11 
and 0.21 for a = 3 and 5, respectively.). This is an expected result when we compare the path- 
loss model dependent constants appearing in Theorem |2] which are also numerically computed 
for various values of a for both path-loss models in Table [I] for a comparative illustration. When 
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TABLE I 

/ r°° (t)tdt 

Path-loss Dependent Gaussian Approximation Constant — — T 

W/~G=(t)tdt)2 



Path-loss model Path-loss exponent (a) 

a = 3 a = 4 a = 5 

Gi(t) = 1 , a 1.564 2.3838 3.1688 

(1 + *) 

G 2 (t) = — — 1.0501 1.1972 1.2713 



we compare the effect of different path-loss models on the Gaussian approximation, we observe 
that the match between the simulated WMAI distributions and the normal distribution is slightly 
better for G 2 (t) (e.g., when A = 0.1 and a = 5 without fading in Fig. [2l the maximum deviation 
between the simulated WMAI distribution and the normal distribution is 0.3 and 0.21 for Gi(t) 
and G 2 (t), respectively.). This is also an expected result when we compare the path-loss model 
dependent constants appearing in Theorem |2] (i.e., see Table H again). 

V. Performance Bounds: Outage Capacity and Sum Capacity 

In this section of the paper, we present further applications of our Gaussian approxima- 
tion bounds above to derive various performance limits and metrics for wireless networks. In 
particular, we will illustrate two potential applications of our results to bound the single link 
outage capacity and the ergodic sum capacity for spatial wireless networks in the current and 
next subsections. These capacity measures can be computed in an exact form for Rayleigh 
fading channels [22] but there does not exist such closed form expressions for general fading 
channel models. Similar to our analysis above, we will employ two different path-loss models 
G\{t) = (i+ t \ a an d G 2 (t) = with various values of a. Similar conclusions continue to hold 
for other path-loss models. 

A. Single Link Outage Capacity 

We start our analysis with outage capacity calculations. To this end, we introduce more 
structure and modeling parameters into the communication set-up under consideration to express 
the problem of bounding the single link outage capacity as general as possible. 
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1) Problem Set-up: We consider a planar wireless network that contains a test transmitter- 
receiver (TX-RX) pair whose separation is d [unit distance] (usually in kilometers). Without 
loss of generality, the test receiver node is assumed to be located at the origin 0. The wireless 
communication link between the TX-RX pair is subject to fading, path-loss, background noise, 
and interference signals emanating from other transmitters in the network. Hence, the maximum 
(short-term) data rate of reliable communication supported by the victim link, which is to be 
defined shortly below, fluctuates as a function of all such effects. We assume that the TX encodes 
data at a rate R nats per second per hertz, and the outage event occurs whenever the maximum 
rate supported by the victim link is below R. This formulation is appropriate for delay sensitive 
traffic for which the communication delay requirement is smaller than the time scale of channel 
variations. 

Let iVo be the power of the (complex) Gaussian background noise present at the RX. We 
define the signal-to-noise-ratio (SNR) of the communication system under consideration as the 
ratio SNR = where P, as defined in Section UH is the transmission power common to all 
transmitters. With a slight change of notation, we will denote the total interference power at 
the RX as I\(P) to simplify the further notation below. Further, we will denote the cumulative 
distribution function of I\(P) by F\. The same assumptions in Section [EI] continue to hold for 
the spatial distribution of interfering transmitters here. 

Then, the signal-to-interference-plus-noise ratio (SINR) at the RX is given as 

SINR = f G( ? - - , (9) 

SNR + p^tJa(I) 

where PG > 1 is the processing gain of the system, and H is the fading coefficient, with 
a finite mean value, for the link between the TX-RX pair. Roughly speaking, PG = 1 case 
represents a narrowband communication scenario, whereas PG > 1 case signifies a broadband 
(e.g., a CDMA network) setting. For a given realization of channel states, the maximum rate 
of reliable communication supported by the victim link in nats per second per hertz is equal 
to log (1 + SINR), and therefore the communication between the TX-RX pair is said to be in 
outage if log (1 + SINR) < R. For a given target outage probability 7 G [0,1], the outage 
capacity for the victim link is defined as 

C A ,outag C (7) = sup {R > : Pr {log (1 + SINR) < R} < 7} , (10) 
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which is the maximum data rate supported by the victim link with outage probability not 
exceeding 7. 

2) Theoretical and Simulation Results: The next theorem provides the upper and lower bounds 

On CVoutage (7). 



Theorem 3: Let £ (h, R) 



^f-SNR-MPG-E[J A (l)] 



VVar[/ A (l)] 
C x . outage 

(7) is upper and lower bounded as 



CA,outago (7) > sup \ R > : 1 - E 
and 



. For the communication scenario above, 



max \ 0, Q AjOUtagc [C[H,R))\lf r ^ KXR 



|^ G(d) 



l(e«-l) \ 



C 



A, outage 



(7) < sup { R > : 1 - E 



min < 1 , Q 



'A.n„,,„> \C [H, R) ) f LJ„ .SNi 



G(d) 



<7MH) 



< 7 ^ ,(12) 



where Q } 



fy(x) ± and and c(x) are as given in Theorem [TJ 



'A,outageV A 7 ^ - 1 - ' 

Proof: Let g be the probability density function for H, possibly a different density function 
than q. By straightforward manipulations, we can express the outage probability Pr (Outage) = 
Pr{log(l + SINR) < R} as 

Pr (Outage) = Pr {SINR < e R - l} 

= ^ 00 pr|(^[-SNR-^PG</A(l)}g(/ i )d/ i 

PrW^ 



1 - 



SNR- 1 (c- R -l) 



SNR PG ^ q(h)dh 



where the last equality follows from the fact that is a positive random variable, and we 

have (^1 - SNR- 1 ) PG < if and only if h < SNR ^ff ~^ . By using Theorem [TJ and the 
natural bounds and 1 on the probability, we can upper and lower bound Pr (Outage) as 



Pr (Outage) < 1 



.^-i^-Q ^ {°> ^.outage (C (K R))} 
GP1 



1 - E 



max \ 0, g AiOUtagc (C [H, R) ) \ l/--. > sx 1 



«('() 



and 



Pr (Outage) > 1 



SNR- 1 (e iJ -l) 
G{d) 



min {l,Q\~ 

outage 

(((h,R))}q(h)dh 



1 - E 



min < 1 , Q 



,outage 



C(H,R 



L f ~ 3NR-1 (,«-!) 
\ - G(d) 
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Fig. 3. Upper and lower bounds on Ca. outage (7) for the path-loss functions Gi(t) 

(bottom figures). Transmitters are distributed over 

effects are modeled by using the Nakagami-m fading model with rn parameter set to 5. For this model of spatial distribution 

of interfering transmitters, c(x) is give: 

a = 4, d = 1, PG = 100 and rj = 0.5.) 



(l+i)° v~r "o "i — i +t o 

2 — B (0, rjd) according to a stationary PPP of intensity A. All fading 



of interfering transmitters, c(x) is given as c(x) = -L- mg3 3 ^ c3(t)tdt m i n fo.4785, Ixril") ■ (SNR = 20 [dB], 



("• H a) J (J^G=(t)tdt)2 



The proof is completed by observing that the upper (lower) bound on the outage probability 
crosses the target outage probability 7 earlier (later) than Pr (Outage) as R increases. ■ 
We plot CA )0 utage (7) and the corresponding bounds in Theorem |3] as a function of A in Fig. 
[3] for both path-loss functions. For this simulation study, we assume that interfering transmitters 
are uniformly distributed over IR 2 — B (0, rjd) according to a stationary PPP of intensity A [nodes 
per unit area], where rj G [0, 00). Here, B (0, rjd) can be interpreted as the exclusion zone around 
the RX in which no interfering transmitters are allowed, perhaps due to physical limitations 
(i.e., small values of rj), or the operation of the MAC layer (i.e., CSMA-CA protocol), or the 
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interference cancellation property [|23l . 7/ is just a modeling parameter that allows us to control 
the radius of the exclusion zone, and the familiar stationary PPP model without any such holes 
can be recovered by setting it to zero. 

We set SNR to 20 [dB], and the path-loss exponent a to 4. All fading effects are modeled by 
means of the Nakagami-m fading with unit mean power gain and m — 5. Similar observations 
continue to hold for other parameter selections. We observe that our bounds closely approximate 
CA,outago (7) within one nats per second per hertz for moderate to high values of A. Especially, 
for high values of A, our bounds become very tight, and they almost coincide with C> )OUta ge (7). 
Considering the simulation and analytical results presented in Subsection [TV] this is an expected 
result since our Gaussian approximation for the interference power at the RX becomes more 
accurate in the dense network limit. 

3) Scaling Behavior: Our results could also provide a high level perspective about the detri- 
mental effects of the network interference on the single link outage capacity. With a slight abuse 
of notation, let 7 (h, R) = 1 -F x (pG 0^ - SNR" 1 ) V 7 (h, R) is the conditional probability 
of outage given H = h when the TX encodes data at rate R. For simplicity, we will assume that 
the TX can track variations in H, and adjust the rate of communication R\ (h) as a function of 
observed values h of H . For a given target outage probability 7 > 0, a reasonable rate se 
policy is to choose R\ (h) such that the outage probability is always equal to 7 for all h 
for any e > 0, we can find A such that 

/PG y^-SNR-^-E^l)]' 
l-7-e<* ^ ' I < 1-7 + e 

for all h > and A > A since the interference distribution can be uniformly approximated by the 
normal distribution. This implies log I 1 H — —, hG \ d ^ r I < R\(h) < 

b \ SNR- 1 + I l J (E[/ A (l)] + v /Var[/ A (l)]vl>-i(l-7+e)) J ~ K ' ~ 

log f 1 H hG ^ = r I , uniformly in h for all A large enough. Hence, 

6 V SNR- 1 + F |;(E[/ A (l)] + ^Var[7 A (l)]^-i(l- 7 -e)) J 3 5 5 

Rx (h 



ection 
Then, 



we conclude that E R\ \^Hj = (^) as A — > oo. Similar but slightly more involved arguments, 
which are sketched in Appendix El also show that C\outagc (7) = © (x) as A — > 00. Indeed, 
this is exactly the behavior observed in Fig. |3] For example, a fivefold increase in A from 20 to 

'Optimum rate selection policy maximizing E (^)j subject to E ^7 (^H,R\ (^))] < 7 is an infinite dimensional, 
possibly non-convex depending on F\, functional optimization problem whose solution is out of the scope of this paper. 
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100 nodes per unit area results in a fivefold decrease in C^outage (7) from 0.25 to 0.05 nats per 
second per hertz in the bottom righthand side performance figure in Fig. |3] 

B. Sum Capacity for Spatial Multiple Access Networks 

We now illustrate another application of our results to bound the capacity of a spatial multiple 
access network in which transmitters are distributed according to a PPP with intensity parameter 
A, and they all transmit to a common base station (BS) located at 0. The same assumptions in 
Section HI continue to hold for the spatial distribution of transmitters here. This set-up requires 
us to interpret the test receiver node above as the common BS, and I\(P) as the useful signal 
power for information flow from transmitters to the BS. Using the same notation above, the 
ergodic sum capacity of the network in nats per second per hertz is equal to 



C A , sum (SNR) 



log ( l + SNR^# fe G(||X 



k\\ 2 ) 



k>l 



E[log(l + / A (SNR))] ) (13) 



which is achievable by using complex Gaussian codebooks and successive interference cancel- 
lation receiver [j24]. Implicit in this formulation is that the communication delay requirement is 
much longer than the time scale of channel variations (i.e., delay insensitive traffic) so that the 
BS can average over the fluctuations in the channel to achieve the communication rates in (fT3l . 
The next theorem provides the upper and lower bounds on C^sum (SNR). 

Theorem 4: For the communication scenario above, C^sum (SNR) is upper and lower bounded 

as 



and 



where Qx( x ) = ^(x) ± ^7=, and \I/(x) and c(x) are as given in Theorem [TJ 
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c(x) is given as c(x) 



7 ff3 Jo° G 3 (t)tdt 

H 2)i (/ °°G=(t)tdt)i 



min 0.4785, 



(SNR = [dB] and a = 4) 



Proof: Since log (1 + I\ (SNR)) is a positive random variable, C A sum (SNR) is equal to 

roc 

CA,s Um (SNR) = / Pr{log(l + / A (SNR)) > a;} da; 
Jo 

= Pr{/ A (SNR) >e x -l}dx 

! Pr ( J a (SNR) - E [I X (SNR)] < g - 1 - E [/ A (SNR)] ] ^ 
\ VVar [I x (SNR)] " ^Var [J A (SNR)] J 

We complete the proof by substituting the bounds in Theorem [T] in (fl?)) and using the natural 

bounds and 1 on the probability. ■ 
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We plot C\sum (SNR) and the corresponding bounds in Theorem @]as a function of A in Fig. [4] 
for both path-loss functions. For this simulation study, we assume that transmitters are uniformly 
distributed over M 2 according to a stationary PPP of intensity A [nodes per unit area]. We set 
SNR to [dB], and the path-loss exponent a to 4. Fading effects are modeled by means of 
the Nakagami-m fading with unit mean power gain and m — 5. Similar observations continue 
to hold for other parameter selections. Again, we observe that our bounds closely approximate 
CA,sum (SNR) within one nats per second per hertz for moderate to high values of A. Especially, 
for high values of A, our bounds become very tight, and they almost coincide with C\ ^ m (SNR). 
With increasing values of A, it is observed that CA, sum (SNR) grows logarithmicallyl which is 
also a property demonstrated by the derived bounds. 

VI. Conclusions 

In this paper, we have focused on the statistical characterization of wireless multi-access 
interference distributions when transmitters are spatially distributed over the plane according to 
a Poisson point process, which is not necessarily stationary. The signal propagation at the physical 
layer is modeled by means of a general bounded and power-law decaying path-loss function. 
Other random wireless channel dynamics such as fading and shadowing are also incorporated 
in the employed signal propagation model. 

We have shown that the wireless multi-access interference distribution converges to the Gaus- 
sian distribution with the same mean and variance at a rate where A is a modeling parameter 
controlling the "intensity" of the planar Poisson point process generating transmitter locations, 
and c(x) > is a function which depends on the shape of the path-loss function and on the 
point x G R at which we want to estimate the interference distribution. c(x) decays to zero at 
a rate |x|~ 3 as the absolute value of x increases. This behavior makes the derived bounds tight 
for any finite value of A. An explicit expression for c(x) appearing in our approximation results 
has also been provided, i.e., see Theorem CD 

We have performed an extensive numerical and simulation study to verify the derived theo- 
retical bounds. A very good statistical match between the simulated (centered and normalized) 
wireless multi-access interference distributions and the normal distribution with zero mean and 



4 Expected total received power at the BS grows linearly with A, i.e., see j I St in Appendix Icl 
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variance one has been observed even for moderately small values of A. Since there are no closed 
form expressions available for the wireless multi-access interference distributions under general 
bounded path-loss models at the present, these results are expected to help researchers in the 
field by simplifying the derivation of closed form expressions for various performance bounds 
and metrics for wireless networks. In particular, applications of our Gaussian approximation 
results have been illustrated to derive tight upper and lower bounds on the outage capacity for 
a given victim link in a spatial wireless network and those on the sum capacity for spatial 
wireless multiple-access networks. It has been observed that the derived performance bounds 
can approximate these capacity metrics up to one nats per second per hertz for moderate to high 
values of A. With increasing values of A, the approximation precision improves significantly, 
and the derived performance bounds almost overlap with the simulated outage capacity and sum 
capacity curves. 
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Appendix A 
Proof of Lemma CD 

We can find constant B 1 > and K > such that G(t) < Kt~ a for all t > B x . Then, 

(1 _ e -'PhG(t)j q (h) p (t)dtdh 

poo f*B\ poo poo 

< / (1 - e - sPhG ®) q(h)p(t)dtdh + / / [1 - e- sPKht ~ a ) q(h)p(t)dtdh 

JO JO Jo Jbx ^ ' 

/B\ poo poo 

p(t)dt + J J (l-e- sPKht ~ a )q(h)p(t)dtdh. (15) 

The first integral in ([TBI is finite since A is locally finite. To show the finiteness of the second 
integral, we divide it into two parts as follows. 

J" (l - e - sPKht - a ^j q(h)p(t)dtdh 

p p K POO POO POO 

= I / f 1 - e - sPKht ~ a j q(h)p(t)dtdh + / f 1 - e"™") q(h)p(t)dtdh. (16) 

^0 J Bi J „„ J B\ 



JO 
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The first integral in ([TBI can be bounded as 

/ aPK / (l - e - s ™" Q ) q(h)p(t)dtdh < (l- e -*~°) p(t)dt, 

which is finite since 1 — e~* = O (t~ a ) and pit) = O (i a ~ 1-e ) as t — > oo. Hence, proving the 
finiteness of f°°i ( 1 — e - sPKht ~ a \ q[h)p{t)dtdh will complete the proof. To this end, we 

sPK 1 \ / 

need the following lemma. 

Lemma 5: 1 — e~ at < 2a (1 — e~ a ) t~ a for all a > 1 and t large enough. 

Proof: We let / t (a) = 1 — e~ a * " and <7t(a) = 2a (1 — e~ a )t~ a . For a = 1, we have 
lim^oo t^j- = 1 and lim^oo = 2 (1 — e _1 ) > 1. Hence, there exists a constant £? 2 > 
such that gt(l) > for all t > B 2 . Fix an arbitrary t greater than B 2 . Then, 

= t-«e-«t-° and ^ = 2T« (l + ae"° - e"») . 
da da 

Thus, <7t(a) grows faster than f t (a), implying that g t (a) > f t (a) for all a > 1 and t > B 2 . ■ 
By using Lemma [51 we can upper bound the second integral in (fT6l) as 

/"°° /°° (l - e -^t-^ q ^) p (t)dtdh 

(*Bs poo poo 

< / p{t)dt+ / / 2sPKh(l-e- sPKh )q(h)r a p(t)dhdt (17) 

for some positive constant B 3 large enough. The first integral in (TTvT > is finite due to local 
finiteness of A. The second integral in (fTTT) can be upper bounded by IsPKmn t~ a p(t)dt, 
which is finite since m H < oo and p(t) = O as t — > oo. 

Appendix B 
Proof of Lemma [2] 

We will show that C In (s) converges to C Ix (s) pointwise as n tends to infinity. Since U liTl , . . . , U^^ n 
are independent, we have 

C In (s) = (l - £ jH £ (1 - e~ sPhG ^) p(t)q(h)dtdh) 

We have J °° £ (l - e - sPhG ^) p{t)q{h)dtdh converging to J °° J °° (l - e - sPhG ^) p{t)q{h)dtdh 
as n tends to infinity, and J °° J °° (l — e - sPhG ( t )) p(t)q(h)dtdh < oo by Lemma [TJ Hence, by 
observing that A (IR 2 ) = oo, we have 

(POO POO \ 
-X / (1 - e- sPhG{t) ) q(h)p(t)dtdh\ 
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Appendix C 
Proof of Lemma [3] 

By using Campbell's Theorem lfT8l . we have 



POO 

E [J A ] = \Pm H / G(t)p(t)dt (18) 
Jo 



and 



poo 

Var = \P 2 mjj2 / G 2 (t)p(t)dt. (19) 

Note that our assumptions on G(t) and p(i) ensure that E and Var are both finite. Let 
Uk y n be defined as in the proof of Lemma |2] Let also m fe n = E \PH k G (C4 )Tl )] an d cr 2 n = 
Var [PH k G Then, E [J n ] = |~A n ] mi,„ and Var [J n ] = |~A n ] a\ . We can explicitly write 

mi„ as 







which implies that lim n _^oo E [/„] = E since A (M 2 ) = oo. Similarly, we have 



-l,n~ A 

which implies that lim n ^oo Var [I n ] = Var 

Appendix D 

Normal Approximation Results for a Non-stationary PPP 

Our normal approximation bounds given in Theorem Q] are valid for both stationary and non- 
stationary PPPs. In this appendix, we will illustrate the validity and utility of these bounds for 
a non- stationary PPP To this end, we place the test receiver node to the origin, i.e., X Q = 0, 
and consider a non-stationary PPP, still denoted by $a, with mean measure density f given as 

if II l| 2 > r 



f{x) = { " , (20) 

if \\x\\ 2 < r 

where r > is a given positive constant. As in Section |V] r can be interpreted as the radius 
of an exclusion zone in which no interfering transmitters are allowed, perhaps due to physical 
limitations, or the operation of the MAC layer {i.e., CSMA-CA protocol), or the interference 
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cancellation property 11231 . Since the test receiver node is located at the origin, we have T (x) = 
\\x\\ 2 , and 

A or -1 ([(),*]) = A {x e M 2 : \\x\\ 2 < t} 

71 |j2 l{r<\\x\\ 2 <t}dx 

\\ X \\2 

= 27rA(log(0-log(r)) 

for t > r. Hence, the transformed PPP J2k>i^T(x k ) has the mean measure density p\(t) = 
A^l{ t > r y. Finally, we have the following theorem approximating the WMAI distribution in this 
case. 

Theorem 5: Assume $a is a PPP with mean measure density / given as in (l20l) . Then, for 
all ieK, 



p. h E[/ A ] <x y v] - • 



< ^ (21) 
7 Var[/ A l ~ | ~ VA' 

where V(x) = -4- P e~^dt and c(x) = -± — mg3 3 r^ 3 ^ min ( .4785, M 

Proof: Directly follows from Theorem \T\ after substituting ^-ljt^r} f° r p(t)- ■ 
In Figs. \5\ and [6l we plot our numerically computed Gaussian approximation bounds and 
simulation results, respectively. Since the key messages conveyed by these figures are similar 
to those explained in Section [TV] we do not repeat them here again. However, several remarks 
are in order. Gaussian approximation bounds and simulation results given for Gi(t) resemble to 
those given for G 2 (t) much more closely when compared to such bounds and simulation results 



2 



given in Figs. \T\ and [2l This is because the path-loss dependent constants, i.e., 
% = 1,2, are more similar to each other with this particular choice of non- stationary spatial 
distribution of interfering transmitters over R 2 . For example, this constant is equal to 1.27 for 
Gi(t) and equal to 1.11 for G 2 (t) when a is set to 4. This also explains why we observe a 
significantly better statistical fit between the WMAI distribution and the normal distribution for 
Gi(t) in this case, i.e., see the corresponding constant in Table H Finally, it is also noteworthy 
to mention that the WMAI distributions become closer to the Normal distribution for smaller 
values of r as the path-loss dependent constants decrease when r decreases. 
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Fig. 5. Upper and lower bounds on the centered and normalized WMAI CDFs for the path-loss functions Gi(t) — njtp 
(lefthand side figures) and G2(t) = ^ * „ (righthand side figures). Interferers are distributed over R 2 according to a non- 
stationary PPP with mean measure density given as in d20t . The effect of fading is also illustrated in the bottom figures by 
assuming Nakagami-m fading with m parameter set to 5. (a = 4 and r = 0.5) 



Appendix E 
Scaling Behavior of C7A, outa ge(7) 

In this appendix, we will briefly sketch the proof ideas leading to CA, ou ta g e (7) = © (t) as 
A — > 00. Firstly, we observe that 7 (h, R) is a non-increasing and continuous function of h that 



is equal to one for h £ 



0. 



SNR- 



G(d) 



and approaching to zero as h grows to infinity for 
any fixed value of R under our assumptions in Section UH Hence, we can find h\ £ (0, 00) 
such that 7 (h* x , C\ :OUta , gc (7)) is equal to 7. Secondly, we show that < h — liminf^oo h* x < 
limsup A _ >oc h\ = h < 00. The proof of this assertion is based on proof by contradiction. For 
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Fig. 6. Comparison of the simulated centered and normalized WMAI CDFs with the normal CDF for the path-loss functions 
Ci(^) = (i+t)° (lefthand side figures) and G?2(i) = yqr^ (righthand side figures). Interferers are distributed over R 2 according 
to a non- stationary PPP with mean measure density given as in J20b . The effect of fading is also illustrated in the bottom figures 
by assuming Nakagami-m fading with m parameter set to 1. (r = 0.5) 



example, if lim sup^^ h\ = oo, then it follows that limsup^^ 7 (h, C\ )OUt age (7)) = 7 for 
almost all h with respect to the distribution of H. By using Theorem \T\ this result implies that 
for any given e > 0, there exist positive constants B\ and B 2 such that 



SNR" 1 + ^ (E [h{l)\ + V /Var[/ A (l)]vl/- 1 (1 - 7 + e) 

h* x <B 1 -) f (22) 

SNR" 1 + (E [J a (1)] + VVar^lp- 1 (1 - 7 - e) 



for all A > B 2 . Hence, h < B\ < 00, which is a contradiction. A similar contradiction shows h > 
0. Finally, it follows that log ( 1 + — w G( ?~ r ) < CAouta E e (7) < 

y B \ SNR- 1 +p^(E[/ A (l)] + v /Var[/ A (l)]*-i(l- 7 +e)) J ~ K !J - 

log I 1 H ; — —, G , r ) for all A large enough. 

6 \ SNR-i+pt, (E[/ A (l)] + v /Var[/ A (l)]*-i(l-7- e )) J & & 
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